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Abstract 



In this paper we give a precise mathematical formulation of the relation between Bose 
condensation and long cycles and prove its validity for the perturbed mean field model 
of a Bose gas. We decompose the total density p = /9 s hort + Piong into the number 
density of particles belonging to cycles of finite length (p s hort) an d to infinitely long 
cycles (piong) m the thermodynamic limit. For this model we prove that when there is 
Bose condensation, pi on g is different from zero and identical to the condensate density. 
This is achieved through an application of the theory of large deviations. We discuss 
the possible equivalence of p\ ong ^ with off-diagonal long range order and winding 
paths that occur in the path integral representation of the Bose gas. 
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1 Introduction 



In 1953, Feynman analyzed the partition function of an interacting Bose gas in terms of the 
statistical distribution of permutation cycles of particles and emphasized the roles of long 
cycles at the transition point pQ . Then Penrose and Onsager, pursuing Feynman's arguments, 
observe that there should be Bose condensation when the fraction of the total particle number 
belonging to long cycles is strictly positive These ideas are now commonly accepted and 
also discussed in various contexts in systems showing analogous phase transitions such as 
percolation, gelation and polymerization (see e.g. [I], [H])- However, to our knowledge, 
there had not appeared a precise mathematical and quantitative formulation of the relation 
between Bose condensate and long cycles until the the recent work of Siito jHj dealing with 
the free and mean field gas. The purpose of this paper is to formulate such a relation and 
to prove its validity in a model of a Bose gas having some genuine short range interactions 
between modes, the perturbed mean field model studied in 

In Section 2 we recall the general framework for cycle statistics independently of any specific 
model. The arguments used below have been well known for a long time in various contexts, 
see e.g. Ginibre in jBJ for a proof of the convergence of virial expansions of quantum gases, 
Cornu [H] in relation to Mayer expansions for quantum Coulomb systems, Ceperley [10. 
for numerical simulations on Helium; a pedagogical account can be found in [IT] . Using 
standard properties of the decomposition of permutations into cycles, the grand-canonical 
sum is converted into a sum on cycle lengths. This makes it possible to decompose the 
total density p = p sn ort + Piong into the number density of particles belonging to cycles of 
finite length (p s hort) and to infinitely long cycles (piong) in the thermodynamic limit. It is 
conjectured that when there is Bose condensation, pi on g is different from zero and identical 
to the condensate density. The main result of the paper is to establish the validity of this 
conjecture in the perturbed mean field model. 

This model is diagonal with respect to the mode occupation numbers with a non mean 
field interaction between them. It has been thoroughly analyzed in jT2j, |E] and it is 
revisited here in Sections 3 to 6 from the view point of cycle statistics. In order to compute 
Piong = P~ Pshort, we find it useful to introduce a modified statistical ensemble equipped with 
an additional external field that gives an extra weight to cycles of finite length. Then p sn0 rt 
is related to the derivative of the pressure with respect to this external field. In Section 
3 the mean perturbed field model is presented, and an easy explicit calculation using the 
modified ensemble establishes the identification of pi on g with the Bose condensate for the free 
gas. Section 4 is devoted to the formulation of the variational principle that leads to the 
same identification for the interacting gas. The justification of this variational principle is 
provided by an application of the theory of large deviations (Section 5) and the mathematical 
proofs are found in Section 6. In the concluding remarks (Section 7) we discuss the possible 
equivalence of pi on g 7^ with off-diagonal long range order and winding paths that occur in 
the path integral representation of the Bose gas. 



2 Cycle statistics 



2.1 Cycle representation of the partition function 



We consider a system of n identical bosons enclosed in a region A in thermal equilibrium at 
temperature /c#T = The states of these bosons belongs to the symmetrization 7Y™ symm 
of the n-fold tensor product of Hf n of the one particle Hilbert space H A and their total energy 
is given by a symmetric Hamiltonian Ha- The projection -P S y mm of Hf n onto symm reads 

where U n : 7Yf™ h- > Hf n is a unitary representation of the permutation group S n on Tif n 
defined by 

[7^(01 (8) 02 <E>- • = 07r(l)®07r(2)<S>---®07r(n), <frj E H A , j = 1, . . . , 71, 7T G 5 n . (2.2) 

The canonical partition function can be written as 

Z2 = trace ^ _ [e"^] = trace ^ [P^fi^] = 1 ^ ^ ( 2 ' 3 ) 

7TGS„ 

where / a (tt) is a function on defined by 

/ A (?r) = trace w ®n [[^e - ^] . (2.4) 



Each permutation tt £ S n can be decomposed into a number of cyclic permutations of 
lengths qi,q2, ■ ■ ■ ,q r with r <n and q± + q 2 + ■ ■ ■ + q r = n. We consider the set Q = U r& ^fl r 
of unordered r-tuples of natural numbers q = [qi, q 2 , ■ ■ ■ , q r ] € Q r for r = 0, 1, 2, . . ., and let 

— 1\ + ?2 H h g r for q G 11. Then a decomposition of n E S n into cycles is labelled by 

q G f2 with |q| — n. We recall the following facts on the permutation group. 



• The decomposition into cycles leads to a partition of S n into equivalence classes of 
permutations C q , |q| — n. 

• Two permutations n' and n" belong to the same class if and only if they are conjugate 
in S n , i.e. if there exists ir e S n such that 

n" = tt-Vtt . (2.5) 



• The number of permutations belonging to the class C q is 



ni 



Wq!(5l52 • • - 5r) 

with n q ! = rii\n 2 \ • • • rij\ ■ ■ ■ and rij is the number of cycles of length j in q. 



(2.6) 



We observe that for a symmetric Hamiltonian ([H A , = 0, 7r G /a(tt) is constant on 
C q . Indeed for tt', n" G C q , one has by (12. 5|) 

/ a (tt") = trace w ®n [^»e~^ A ] = trace H ®n [U^U^U^e'^] 

= trace w ®« [[/^[/^"^[/J = trace w ®« [t^e - ^] = / a (tt') . (2.7) 

Therefore / A (7r) only depends on the cycle decomposition q of tt, and we shall write from 
now on 

/a(tt) = / A (q), tt G C q . (2.8) 

For instance / A (q) = A^q) can be evaluated on the permutation 7r q = (1, . . . , qi)(qi + 
1, . . . , q% + 52) . . . (n — q r , . . . , n) representative of C q . Let U q be the operator representing 
the cyclic permutation (1,2, ... ,q). Then 

U^ = U qi U q2 ---U qr , (2.9) 



/ A (q) = trace nfn [U gi U q2 ■ ■ ■ U qr e~^] (2.10) 

and using (|2.6|) 

_ — . Tt I 

E AW = ^7— \ A(q) • (2.11) 



TTSCq 



Then, summing first in (j2.3|) on the permutations belonging to a class C q and then on all 
possible classes gives 



^E E E« 



rvn 

f=l qSf2 r , |q|=n ir£C q 

= hth E v Eaw 

r=l q£N r ,|q|=n 7reC q 
n 1 1 

= E^ E 7 <&,■■■,<&•) (2.12) 

In the second line we have taken into account that replacing the summation on unordered 
r-tuples in Q r by ordered r-tuples in N r , we are over counting by a factor r!/n q !. Finally 
to get rid of the constraint of fixed particle number n, we introduce the grand-canonical 
partition function with chemical potential [i 

00 00 r n m 

n=l r=l ' qgN r i=l ^ 

This leads to define the grand-canonical distribution of cycles on the space Vl — U f .gj^r2 r of 
r-tuples of arbitrary integer length by the probability measure 

1 1 ' e Pm 

P£(q) = J] Q2,-..,q r ), qeH, (2.14) 

OO 

E E P^(gi,?2,...,gr) = i 



One then introduces the cycle density and cycle correlations in the usual way. Consider the 

1 r 

observable — 5 qSi i.e. the number of cycles of length q per unit volume. Its grand- 

I I i=i 

canonical average with respect to P A is given by 

\ ' ' i=i i 



= \a\z£— ^^~Ty. 2. 11— /A(ff,ft.-.ff(r-i))- (2-15) 

The quantities 

1 °° 1 <9 

P,M = M h lS ; and = = __l»aj ( 2 . 16 ) 

are the usual grand-canonical pressure and particle number density of the gas. 



2.2 Particles in long and short cycles 

We first note that expressing the partition function in terms of cycles corresponds to a 
partitioning of the particles into groups of the same size as the cycle lengths. It is natural 
therefore to ask for the probability p A (q) for a particle to belong to a cycle of length q (we 
drop the chemical potential p from the notation and it is understood that p and the inverse 
temperature j3 are kept fixed). Since the average number of particles to be found in g-cycles 
is |A|gp A (g) and the total average particle number is |A|p A this probability is 

VM = f> A (g) = 1. (2.17) 

q =i 

Let us assume at this point that infinite volume limits of the cycle and total densities exist, 
namely 

lim p A (q) = p(q), lim p A = p . (2.18) 

A— >oo A^oo 

Then liniA^ooPA^) = Qp(q)/p — p(o) is the probability to find a particle in a g-cycle, 
q < oo, in the infinitely extended system. Hence the fraction of particles pi on g/p belonging 
to infinitely long cycles is 



l-$>(ff) (2-19) 



p 

The conjecture [2| is that pi on g is different from zero when there is Bose condensation and 
identical to the condensate density. One also says in this case that the probability p(q) is 
defective or that cycle percolation occurs. 

From the above definitions one can write 

oo oo oo 

Plong = P~Y] Wig) = lim V qp A {q) - V qp{q) 

' ' A — >oo — — 

q=l q=l q=l 

oo oo 

= lim ?Pa(?)~ J2 9P(0) (2-20) 

A->oo z — ' — ' 
q=Q+l q=Q+l 



for any Q > 0. Thus letting Q — > oo and taking into account that the last sum in (|2.2()jl is 
convergent one has 



pi ong = lim lim V qp A (q) (2.21) 

Q~ >oo A-^oo ' * 
q=Q+l 



showing that pi on g is carried by the tail of the cycle number density as A — > oo. 
One can alternatively consider the quantity 



Pshort = P - Piong = hm lim y^qp A (q) (2.22) 



■»00 A— KX) 

9=1 

which is expected to coincide with the density of particles that do not belong to the conden- 
sate and can be studied more easily. It will turn out to be convenient to define a modified 
partition function depending on a parameter A > by 

S q ,a,a = 1 + ^- II : /a(?i, &,■■■,?»■) (2.23) 



rGN ' q€N r i=l 



where 



and the corresponding modified pressure 



Pq, a, a — hi S Q) a, a- (2.25) 



One immediately sees that 

E^) = ;^W A=0 - (2.26) 



9 



<9A 

q=l 



Thus if P Qj a = hmA^oo P Qi a, a has a thermodynamic limit, the convexity of P Qi a, a(p) implies 
the convergence of the derivative jHj and 



d 

p short = lim lim \^qp A (q) = lim ^t-P q ,a • (2.27) 

Q^oo A^oo Q^oo OA A=0 



Note that the order of limits is crucial here. The procedure is the same as generating the 
density from the pressure by functional differentiation with respect to an external field. Here 
the "external field" X0q(q) gives an additional weight to cycles of length q < Q. 



3 The perturbed mean field model 

The bosons have mass M and are enclosed in a cube A C M. d centred at the origin with 
periodic boundary conditions. The one-particle space H A = L 2 (A) is generated by the basis 
of plane waves with wave number k E A* = {27r//|A| 1//d | / G Z d } where |A| is the volume of 
A. In second quantized form, the Hamiltonian of the perturbed mean field model expressed 



in terms of the boson creation and annihilation operators a* k and a k , k G A*, [a k , a* k ,} = 5k,k'- 
reads 

H K = T K + -£-N\ + ±- J2 v(k,k')N k N k , (3.1) 



2V IV 



where N k = a* k a k are the occupation number operators of the modes k, and 

T A =J2 < k ) N k, <k) = h 2 \\kf/2M, N A = J2 N k (3-2) 

fceA* fceA* 

are the kinetic and total particle number operators. It differs from the standard mean field 
model by the last term of ()3.1|) which introduces the coupling kernel v(k,k') between the 
modes. We assume that v(k,k') is a bounded, continuous and positive definite kernel with 
suitable decay properties. 

Since the perturbed mean field model is defined in terms of the mutually commuting occupa- 
tion numbers N k , it can be studied with the techniques of classical probability theory. The 
iVys form a set of classical random variables and the equilibrium properties are expressed in 
terms of their distribution functions. In particular, in the thermodynamic limit, the average 
occupation number distribution is a measure v(dk) which is continuous except for a possible 
Dirac 5 at k = corresponding to a macroscopic occupation of the k = state. Thus when 
there is a non vanishing condensate density u c , the equilibrium distribution is given by 

dk 

v(dk) = u c 5(k)dk + Ve{k)j^y d (3.3) 

where the continuous part v e (k) describes the density of particles in excited states. The 
main result of the paper consists in establishing that the average number density of short 
cycles (J2.22|) is identical with the total density v e = f v e (k)-[0^ of exited particles, and 
consequently that pi ong (I2.21J) indeed coincides with the condensate density of particles v c in 
the k = mode as given in ()3.3j) . 

For a number of specific interaction kernels it is known that Bose condensation occurs for 
a range of chemical potentials. In particular in the case of a Gaussian kernel v(k,k') = 
t! e~ c ll fc-fc II there is a value /i such that if [i G (0,yU ), then p c is non zero [12] . On the 
other hand for the kernel v (k, k') = t>oe _c " fc_fc " there is condensation for all positive \x [T3*] . 
Our present results therefore imply that in these cases the Bose condensate is carried by 
particles in long cycles. 

We compute the function f A (q 1 ,q 2 , . . . , q r ) for the perturbed mean field model ()3.1|) by cal- 
culating the trace ()2.4j) in the plane wave basis 

AM = (^l)---)^n| e ^ A |^7T (1) 5 • • • ; ^7T (n) ) 

ke(A*) r 

{/ n 1™ \l n 

-p e ^ + ^A\ n2 + ^tat S k ^ r II 

\j=l ' ■ 11 j,j'=l / J j=l 



(3.4) 



giving 

A 92, • • • ,q r ) 



ke(A*) r I \i=l 



(3.5) 



As a preliminary simple calculation, we consider the modified partition function S° A A of 
the free gas. In the absence of interactions the function (|3.5J) reduces to 

/a(?i, Q2, ■ ■ ■ , q r ) = E exp \ -P^ £ (h)<li > • (3.6) 
ke(A*)<- I i=i J 

Inserting this into ()2.23|) we obtain the modified partition of the free gas in the form 

-° q ,a, a = i+£^ E n^— — E «J-/?2>(**)» . (3-7) 

reN ' qGN r i=l Vi ke(A*) r I i=l J 

Here the chemical potential, denoted by a, has to be negative, and the corresponding pressure 
is 

P q,a,a = 1iiH q,a,a (3-8) 

This partition function is readily calculated, permuting product with the q and k summations 
one obtains 



!+ E^ E En 

reN kg(A*) r q€N r i=l 



^0 , \ - ' \ - TT exp{/5(o!- efc) + \0q 

""Q, A, A 



exp 



5> QjA (a- £ (A;))l (3.9) 

fcGA* J 



.fceA 

where for y < and A > we have introduced the function 



9=1 ^ 9=Q+1 ^ 

This yields the modified pressure (12 .25(1 in the infinite volume limit 
Notice that for A = 0, 

7r Q)A=0 (y) = n (y) = p-'Hl - e-?y), (3.12) 

so that 

Pl™ = r 1 ]^Ml-e- mt) - a) ) (3-13) 



reduces to the well known pressure of the free gas at chemical potential a < 0. Moreover, 
one explicitly checks from (j3.10|) and ()3.1 1|) that 



A=0 



dk 



(27r) d e^( £ ( fe )" Q ) - 1 



(3.14) 



with v e the total number of particles in exited states. Bose condensation is characterized in 
the three dimensional gas by the fact that remains finite for a = (the critical density). 
The relation ()3.14|) also remains valid for a = showing that the equality p s hort — holds 



even in case when the condensate density is different from zero 
of [6 for the free gas. 



This yields again the results 



4 The variational principle 

We set up a variational principle to determine the modified pressure ()2.25|) in presence of 
the interaction. Its precise mathematical justification is provided by the application of the 
theory of large deviations presented in the next sections. Here we give the form of the 
grand-potentials per unit volume <p = u — pp — Ts = —P (u is the internal energy, p 
the particle number density and s the entropy) corresponding to the perturbed mean field 
model and to the modified statistical ensemble ()2.23|) . These potentials will be considered 
as a functionals of all possible momentum distributions. The latter are described by positive 
bounded measures m(dp) on IR d . Their minimization will yield the corresponding pressures. 

Consider first the measures of the form 

1 r 

m^){dp) = jj^^2q i 5(k i -p)dp, (q, k) = {(qx, h), (q 2 , k 2 ) . . . , (q r , K)} (4.1) 

which are the sum of Dirac distributions centered at the wave numbers ki. The particle 

number density is J m^^dp) = Y^l=iQi/\M = n /\M> anc ^ X\ m (q: k )(^) ^ s ^ ne num ber of 
particles per unit volume that have wave numbers in the range A. One can write ()3.5|) as 

/a(q) = ex P( _ ^l A l w (q,k)) 
ke(A*) r 

«(q,k) = (e,«i(q,k)) + ^(l,m {qjk) ) 2 + ^(m (qik) ,i;m (qik) ) (4.2) 

with the notation (g,m) = J g(p)m(dp), (m,vm) = J v(p,p')m(dp)m(dp'), and where W( q ,k) 
is the energy density of particles in a cycle and wave number configuration (q, k). In view 
of ()4.2j) . we define the energy density functional for general measures m(dp) by 

u(m) = (e, m) + — (l,m) 2 + -(m,vm) . (4-3) 

Concerning the entropy density functional, we recall its standard form jTHJ 



S0,A = ~ k BTT7 ^ \- mk lllmfc ~ ( mfc + 1 ) ln ( m *: + !)] 

' ' fceA* 



(4.4) 



for a gas of independent Bose particles having average occupation numbers m&, k e A*. In 
the thermodynamic limit, 



s (m) = lim s 0tA (m) = -k B 

A^oo 



dk 

— -,[m e (fc) lnm e (fc) - (m e (fc) + 1) ln(m e (Ar) + 1)] (4.5) 



only involves the absolutely continuous density m e (k) of the infinite volume mode distri- 
bution m(dk) ( for instance it is easily verfied that a macroscopically occupied mode, say 
rik=o ~ n, does not contribute to the limit ()4.5|) ; for the free Bose gas, it is well known 
that at any finite temperature the total entropy is entirely due to the non condensed phase). 
Adding the chemical potential term, the grand-potential functional of the perturbed mean 
field model reads 

<p(m) = u(m) — m) — Ts (m) . (4.6) 

It is proved in [7] for a class of positive definite interaction kernels v(k, k') that the minimiza- 
tion of (p(m) yields indeed the pressure of the perturbed mean field model, the minimizer 
u(dk) is unique and is of the form (jH.Hj) . 

We now want to define the entropy density functional s Qi \ for the statistical ensemble ()2.23|) 
that has the additional "A#q field" . For this we first consider the Legendre transform 

7r*(t) = sup (t s - 7T (s)) = ts (t) - n (y (t)) (4.7) 

y<0 

of the function Tio(s) (JHH2J)- It is given by 

TT*(t)=ts (t)-Myo(t)) (4.8) 
where yo(t) is the unique solution of the equation 

* = ^to = p^ri ( 4 - 9 ) 

and thus 

i-i- t 1 



y (t) = (3- L In {j^y) > *o(*) = -T[-k B (t \nt - (t + 1) \n(t + 1))] . (4.10) 
One therefore sees from (|4.5jl that the entropy can be expressed as 

/dk 
j^y d K(m e (k)) (4.11) 

When A 7^ 0, this motivates the definition of an entropy density functional s Qt \(m) for the 
modified ensemble (depending only on the density of the absolutely continuous part m e (k) 
of the measure m(dk)) by 



Ts QjX (m) = - J — -p7r* jA (m e (A;)) 



(4.12) 



where 

Ki,x(t) = SU P (ty- Tr Q ,\(y)) = ty Qt x(t) - tt Qi x (y Qi x (t)) (4.13) 

y<0 



is the Legendre transform of the function (J3.10)) and y Q .\(t) is the unique solution in y of 
the equation 

* = <,a(w)- ( 4 - 14 ) 

Finally the total functional (p Q ^\(m) (the grand-potential of the modified ensemble) is taken 

as 

<f Q , \{m) = u(m) - n(l,m) - Ts Q> A (m) (4.15) 
where u(m) and s Qt \(m) are defined by (|4.3|) and (|4.12|) respectively. 



The corresponding pressure is thus given by 

Pq,x = -<Pq,x{vq,x) 



(4.16) 



where v Qt \(dk) is the measure that minimizes (p Qt \(m). For A = 0, v Qj \=o(dk) = u(dk) 
reduces to the equilibrium mode distribution ()3.3j) of the perturbed mean field model. The 
A-derivative of the modified pressure is 



d_ 

dX 



A=0 OA 



A=0 



(4.17) 



since derivation with respect to u Qt \ does not contribute at the minimum and only s Q: a 
explicitly depends on A in (p Qt \. Therefore we have from (|4.12|) 



9 P 

dX Q ' x 



A=0 



dk d 



(2ir) d d\ 



A=0 



(4.18) 



with u e (k)dk the absolutely continuous part of u(dk). From ()4.13|) and (|3.10j) 



d_ 

dX 



and therefore 



From (fO)) 



and hence 



d_ 



d_ 

a=o dX 



P, 



Q, A 



K Q ,x{yo{t)) 



dk 



A=0 



Kn,o(yo(t)) 



A=0 



lim tt' {y (*)) = Tr' {yo{t)) = t, 



Q— >oo 



9 

lim ttt-Pq a 

Q^oo C/A 



A=0 



r/A - 
(2^ 



■v e (k). 



(4.19) 
(4.20) 

(4.21) 
(4.22) 



In view of ()2.27|) this establishes our main result, the identification of the total excited mode 
occupation v e — J j0^ v e{k) with the short cycle average number density p s hort f!2.22|) : 



Pshort • 



(4.23) 



5 The large deviation argument 



We shall use the large deviation techniques developed in jjj, [lEj (see also references in [T2] ) 
to justify the variational calculation presented in the previous section. Here we present the 
heuristic argument, the technical details are given in the next section. The large deviation 
theory is a generalization of Laplace's method for evaluating the asymptotics of integrals of 
the type J dp(x)e' y ^ x ^ with dp(x) a probability measure and 7 a large parameter. First we 
have to cast our problem in the form of a Laplace integral. 

Let Q rA to be the set of unordered r-tuples (q, k) = {(qi, fci), (#2, k 2 ) . . . , (q r , K)} in N x A* 
and let Q A = \J re ^Q r>A . We introduce a probability measure on Q A associated with the 
modified ensemble (J3.7)) of the free gas at some negative value a of the chemical potential 



1 1 A e^ +A M<?0) 9l 

K, a, M, k) = ^r— - J n ex P 



i=l 



a<0. (5.1) 



From 1)2.23)1 and (J3.5j) we can now write the total partition function (normalized with respect 
to the free one (|3.7jl ) as an average with respect to P° A A of the chemical potential and 
interaction contributions 

ir^= E ^,A, A (q,k)exp(/3|A|G(m (q , k) )) (5.2) 

"Q, A, A q,k€f2 A 

with 1 ^ 

G(m) = (fi — m) — 2 <3 (-'-' m ) 2 — ^ m " 1 vm ^ (5-3) 

It is understood in the sequel that quantities pertaining to the free gas (with index 0) are 
taken at an arbitrarily fixed negative value a of the chemical potential (since the latter 
quantities are not defined when a is positive). The complete partition S Qi A, A function 
depends on the actual value /1, /1 7^ a, of the chemical potential. This produces the term 
(/i — a)(l,m) in (|5.3)l . The a dependence will drop out in the final calculation. 

We observe that the function to be averaged in (|5.2|) depends on q, k only through the 
measure m( q ,k)- We can define a measure K^ X A (dm) on E, the space of positive bounded 
measures on M. d , through the relation 

/ F(m)K° Q>A)A (dm)= £ P QiAiA ( q ,k)F(m (q)k) ) (5.4) 

so that in terms of this measure the expectation ()5.2|) takes the form of an integral on the 
space E 

= / exp(/3|A|GV))K° Q)AiA (dm) (5.5) 

' 'Q, A, A JE 

Then the pressure of the modified ensemble, using (15.5)1 . can be expressed in the form 

= lim -i-ln [ exp([3\A\G(rn))K,x,A(drn) + P Q>x (5.6) 



where the last term is the free (modified) pressure ()3.11j) . 

In a mathematical sense which will be made precise in the next section, the measure 
K° A A (<im) on E behaves or large A as 

K Q)A)A (dm) « exp (-/3|A|/° A (m)) (dm) (5.7) 

where I® x( m ) * s called the rate function. Introducing this asymptotic behaviour in ()5.5|) . 
we see that the partition function S Qi a, a has the form of a Laplace integral (up to a normal- 
ization). Then by a generalization of Laplace's Theorem (Varadhan's Theorem) we get 

P Q , A = sup (G(m) - 1° A (m)) + P Q ° A . (5.8) 

We still have to determine the rate function J° A . One can guess it with the help of a well 
known argument. For a test function / on M v consider 

C° (/) = lim -Lin f exp(0|AK/,m))K° Adm) (5.9) 



Q > AW ' a^c (3\A\ J E ^ ^ 11 w ' ' Q> 
where (f,m) = J f(k)m(dk). Applying now Varadhan's Theorem to ()5.9|) one obtains 

C° A (/) = sup((/,m)-/° A (m)), (5.10) 

m 

that is C® a is the Legendre-Fenchel transform of J° A . Inverting this transform one gets 

J° A (m) = sup((/ ) m)-^ iA (/)). (5.11) 
/ 

We observe that the integral occurring in ()5.9|) has the same form as that in ()5.5|) with 
G(m) replaced by (/, m). Coming back to expectations written in terms of the probability 
Pq A A (q, k) this amounts to replacing in (|5.2J) the partition function S QjAiA by that of the 
free gas ()3.7j) with shifted energies e{k) — f(k) 

exp (/3|A|(/,m))K° jAjA (rfm) 

1 ~ i J_ e ti{<*+ue Q { qi )) qi f ' 1 

— Eh E II r E exp -/3^(^)-M)) % (.5.12) 

"q, A, A r=1 • qgN r i=1 •« ke(A*) r I i=l J 



Cq,a(/) = / 7^^,A(«-^) + /(fc)) + P c ° A (5.13) 



By the same calculation that led to (j3.11|) one eventually gets 

dk 

To find the inverse Legendre transform of C® x (f) it is useful to first single out the singular 
part m c (dk) of m(dk) leading to (see next section) 

7° A (m) = -/(«- eWK(*) + sup ((/,m e ) - C° A (/)) + P° A (5.14) 

The second term in the r.h.s only involves the density m e (k) of the absolutely continuous 
part of m(dk). With the change g(k) = a — e(k) + f(k) and in view of (|4.13J) it is equal to 

dk f dk 

(a - e(k))m e (k) + sup / — - — - [(g(k)m e (k) - 7r Q ,\(g(k))] 



(a - e(k))m e (k) + j 7r* A (m e (A;)) (5.15) 



Finally, 

I° Q>x (m) = -J(a-e(k))m(dk) + J ^ < A K(*)) + p l x (5-16) 

and from (|5.8jl and (J5.3)) 

/ i i f dk \ 

p q^ = ~ }ni E \(e - fi,m) + -a(l,m} 2 + -(m,vm) + j tt* Qi x (m e (k)) J (5.17) 

is precisely the variational problem for the modified ensemble treated in Section 4. 



6 Proofs 

In this section we establish rigorously the identity (|4.20|) which leads to our main result 
(|4.23J1 . This is done in Theorem 16.11 Recall that for Q e N, y < and A > 0, we have (see 

ED} 

« e /%+A)9 « e (3y q 

9=1 HH q=Q+l ^ 

and 

Q oo 

7T 



Q: 

9=1 g=Q+l 

Here ir' Q x denotes (dir Qt \)/(dy). The function y n' Q x (y) is strictly increasing and its 
range is (0, oo). On the other hand 7r* x (t), the Legendre transform of vr Qj a, is decreasing in 
A since (see fpTTSjl ) 

^{t) = -^(vUt)) = -J2e^ it)+X]q < 0. (6.2) 

9=1 

7r* A is of course convex in t but it is also concave in A as can be seen from the following 

argument: 

We have 

From ()4.14|) we get 

o = ^(v*a(*)) - <A(y*A(t))^|^ (6.4) 



which together with ()6.3|) gives 



/}2,_* ^r(2/Q,AW)< A (!/e,A(t))- f-^(2/o,A« 

^(*) = „ J (f] } ^ '— (6-5) 

<9A 2 7Tq,a(2/q,a(0) 

The right hand side of the last identity can be seen to be negative from the explicit expres- 
sions: 

Q oo 

<x(y) = p^qef^ + p E 

9=1 g=Q+l 

^f(y) = ^(y) = /3E^ A)9 - (6-6) 

9=1 



From Theorems 2 and 3 in ^B] we can deduce the next proposition. We note that condition 
(1.1) in is not necessary; one can replace it with the condition that ir' Q A is invertible on 
(— oo, 0) (see the proof of Theorem 1(a) on page 681). 

Proposition 1 If E is equipped withe narrow topology, then the sequence of probability 
measures K° x A on E (defined in \5.4]) ) satisfies the large deviation principle with rate 
function 7° A (defined in \5.1J$ ) and constants /?|A|. 



Explicitly Proposition 1 means that: 
(LD1) Jq a is lower semi- continuous; 

(LD2) For each b < oo, the level sets {m 6 E : 1° \{m) < b} is compact; 
(LD3) For each closed set C 

lim sup(/3|A|) -1 ln K° x A (C) < - inf /° A (m); 

A-+00 ' m £° 



(LD4) For each open set O 



limsup09|A|)- 1 lBK' 

A^oo 





Q, A, A 



(O) > - ml : J° (m) 



We set inf me 7° \{ m ) = oo. 

From Lemma 4.1 in 7] we know that m \— > G(m) (as defined in (|5.3Jl ) is continuous in the 
narrow topology. Then by Varadhan's Theorem we obtain (|5.17j) : 

P Q>x {fi) = -M€ Q)X {m) (6.7) 
where the functional £ Qt \ : E — > M. is defined by 

£ Q ,\(m) = (e-fi,m) + ^a\\m\\ 2 + ^(m,vm) + J ir* Q X (m e (k)). (6.8) 



We have the following results: 



1- £q,a(0) = 0. 

2. £ Qi a is convex in m and concave, different iable and decreasing in A. The concavity of 
£ Qt a in A follows from the concavity of 7r* A in A. The latter together with ()6.2|) gives 

which implies that 

-\\m\\ < ^^(m) <0. (6.10) 

CM 

since 

o < ^(v*a(*)) < <A(y*A(*)) = (6.ii) 

From the continuity of G Q> a an d the lower semi-continuity of J° A it follows that £ Qi a (m) 
is lower semi-continuous in m. 



3. Suppose that v is symmetric and is continuously differentiable with bounded deriva- 
tives. Then S Q} \ has a unique minimizer u Qt \ whose singular part is concentrated at 
0. The minimizer satisfies the Euler-Lagrange equations: 

e(0)-^ + o|K A || + (wo,A)(0) = 0, (6.12) 
e(k) - n + a\\u Qt x\\ + (vu QtX )(k) = -y Q ,x(v Q ,\,e(k)), Lebesgue-a.e., (6.13) 

where z/ Qi a, e is the density of the absolutely continuous part of v Qj \. This is proved by 
the same arguments as in Theorem 6 in [Zj. Equation ()6.13j) implies that 

v Q ,\,e(k) = n' QjX (-{£(k) - P + a\\i/ Qt x\\ + K,a)W}). (6.14) 



We shall prove the following theorem: 



Theorem 6.1 

dP Q ,x 
d\ a=o 



dk 



^2exp{f3(y (is e (k))q} 



q=l 



dk 

(27T 



id Q 



(6.15) 



Lemma 6.1 For every sequence {X n } of positive numbers tending to zero, v Q ,x n converges 
to v in the narrow topology. Moreover the sequence u Qy x n is tight and v Q ,x n ,e(k) converges 
pointwise to v e {k). 



Proof: First we prove that the set of minimizers {v Qt x | < A < 1} lies inside a compact 
set. To do this we follow the proof of Lemma 5.1 in [Z| using the additional fact that for 
fixed m, £ Qt x{Tn) is decreasing in A. Since £ QjA (0) = 0, £ q ,x{ v q, a) < and thus £ q ,i(^ Q) a) < 
£q,\(v q , a) < for < A < 1 . Therefore 



^,iKa) = ^,iKa) + GK, a ) + P° 1 

1 

2 



< G{v Q> x) < - a)\\v cl x \\ - ^ a \WQ,\\\ 2 + P q,i 



-t^WKaII - {n ~ a)} 2 + —{fi - a) 2 
1 

2a' 



< -{^-af + Pl, (6.16) 



which implies that v Q) x lies in a level set of Ig l which is compact by LD2. 

Let {A n } be a sequence of positive numbers tending to zero and let v* be a limit point of 
{Vq, A n }- Then there is a subsequence A n , such z/ Qi A n; converges to v* . For m G E we have 

Sj.AnjO'Q.An,) < ^A»,W < ^o(m) (6.17) 

and therefore 



Now 



liminf£ QiA (z/ q ,a „.) <5g,o(m). (6.18) 



^Q, A„, (^Q, A„ ; ) — ^Q, 0(^Q, A„, ) - (£q, 0(^Q, A„ ; ) ~ ^Q, A„, (^Q, An,)) 

> ^Q.oK.An,) - A ni ||i/ QjA ||. (6.19) 



and thus since £ Qi o is lower semi-continuous lim inf £ Qj \ n (u Qt \ n ) > £q,o(v*) and so 
£q, o( z/ *) < £q, o( m ) for all m & E. This implies that z/* = z/ Qi0 = v and therefore a„} has 
only one limit point v and thus z/ Qj An converges to v. From ()6. 14|) using ir' Q x (y) < ir' Q (y + X) 
with, vo = sup \v(k, k')\, we have that 

f<J,An,e(fc) < 7Tq,o(-{ £ W -^ + «II z/ Q,aJ| + (w Q , An )(fc) - 1}) 

< 1}) (6.20) 

for some K and all n. For large ||fc||, — fi — VqK — 1 > e(k)/2 and therefore 

^■VeW< e J 2 _ r (6-21) 
Again from (jdl4j) that z/ Qj \ n , e (k) converges pointwise to 

<o(- " ^ + o||i/|| + M(fe)}) (6.22) 

since (vi/ Q> \ n )(k) — > (vu)(k) for each fc. □ 



Proof of Theorem 16.11 

Recall that P Qj a(aO = — £q, a( z/ q, a)- Therefore 

-Pq, a„ (aO - -Pq, o(a0 = £q, o(^q, o) - £q, a„ (^q, a„) 

= (£q,o{Vq,o) -^Q,A n K,o)) + (^Q,A n (^Q,o) ~ £q, A n K, A J ) 

> g Q ,oKo) - g Q ,A w K,o) > -X n d£Q ' X d ^ Q ' 0) | w (6.23) 
In the last inequality we have used the fact that £ Q; a is concave in A. Therefore 



A=0 



(2vr; 



X] eX P{^Q, o(^Q, 0, e(k))q} 



q=l 

Q 



dk 

j^y d Y^eMPVoMkM- (6-24) 



On the other hand 



Pq, \n (aO - p q, o(a0 = £q, oK, o) - £ Q , a„ K, a J 

= (£q, ( Vq, ) - £q, Q, A„ ) ) + (£q, ( ^Q, A„ ) - £q, A n (^Q, A n )) 

< £ q ,oKaJ -^,a„KaJ < -A n 9gQ '^ Q ' Aj | A=A/ (6.25) 

Thus 



T™ p £ -JWf^ exp{ ' 3 (A " + 



(6-26) 



To show that the righthand side of (|fi.2fij) converges to the righthand side of (|fi.l5j) we note 
that 

Q 

^2exp{/3(\ n + y QiXn (v Q ,\ n ,e(k)))q} < Tr' Q ,x n (y Q ,\ n (v Q ,\ n ,e(k))) = v Q ,\ n , e {k). (6.27) 
Therefore by using the tightness of the sequence z/ Qj \ n we can make both the integrals 

dky2exp{(3(\ n + y QiXn (v clXn ,e(k)))q}, / dk Y] exp{/?y Qi (u Qt , e(k))q} 

M 2 >R q =\ A\H 2 >R q =l 

(6.28) 

arbitrarily small by taking R large enough. On || k\\ 2 < R we use the bound 

Q Q 

J2 e MP(K + y Q> XnKxn,e(k)))q} < ^exp{#?}. (6.29) 

g=l 9=1 



for A n < 1 and the Lebesgue Dominated Convergence Theorem. 



□ 



7 Concluding remarks 

The most general way to define the condensate density p c for an interacting Bose system is 
by means of the concept of off-diagonal long range order 

p c = lim lim p A (x,x') (7.1) 

x- x'\ +oc a|— >oo 

where p A (x,x') are the configurational matrix elements of the one-body reduced density 
matrix. For the perturbed mean field model the infinite volume reduced density matrix 
p(x, x') is simply the Fourier transform of the equilibrium distribution u(dk) (|3.3j) so that 

lim p{x,x')= lim I v{dk)e iHx - x>) = v c (7.2) 

\x— x'\— *oo \x— x'\— >oo J 

since the Fourier transform of the smooth part of the mode distribution does not contribute 
asymptotically. We therefore see that the three possible definitions of condensation, namely 
particles in long cycles pi ong , macroscopic occupation of the k = mode v c and off-diagonal 
long range order p c coincide in this model. It may be conjectured that the relation pi ong = p c 
holds under rather general circumstances (e.g. for models with local pair interactions). It is 
therefore of interest to examine the relation of p A (x, x') with the cycle statistics. For this we 
first consider the average of an arbitrary one-particle observable of the form A(n) = ^2™ =1 A4 
with Ai = A a single particle operator acting on the one-particle Hilbert space Ti. A . Its 
grand-canonical average reads 

00 

(A) A = ^Y1 e ^" trace hT [A{n)e-P H ^ n) ] (7.3) 



n=l 



where H A (n) is the n-particle Hamiltonian. We can write as in ()2.4j) 

^ace H1 „ [A{n)e~W] = i £ ffo) with 

/ A A (vr) = trace n ® n [U{u)A(n)e- pH ^ n) ] . (7.4) 



For the same reasons as in (|2.7|) //(i) = /^(q) only depends on the cycle decomposition of 
7r and all the steps leading to (j2.13|) are the same: 

= e h e n V AA(gi ' & ■ ■ ■ > gr) • (7 - 5) 

" A r =l ' qGN r t=l ^ 

With ()2.9|) one can write the function f A (q%, q-z, ■ ■ ■ , q r ) more explicitly as 



/a (9i,92, • • • ,9r) = trace 



(7.6) 



r x trace _,<8 91 



where the notation H A {n) = H A (qi, g 2 , • • • , q r ) indicates that among the n interacting parti- 
cles qj of them belong to the cycle j, j — 1, . . . , r. It is a symmetric function of q%, q^ ■ ■ ■ , q r 
so that when inserted in (|7.5jl and relabelling the dummy summation indices, it contributes 

as 

r x trace H m^ Hmr [U qi U q2 ■ ■ ■ U gr A( qi )e^ H ^'-^] = 

U qi A{ qi )tmce ^...^ (U q2 ■ ■ ■ U qr e^ H ^'-^)\ (7.7) 

Hence we obtain finally (changing r — > r + 1 in the r-sum) 

00 

(A) A = trace H m [A{q)U q R A {q)} (7.8) 

9=1 

where we have introduced the cycle reduced density matrix R A (q) as an operator acting in 
the space Hf q of the the particles belonging to a g-cycle: 

Mq) = -^Eh £ II trace «r ®-®nr ■ ■ ■ Uqr e-^^-^] 

(7.9) 

The reduced density matrix is obtained from (|7.8j) by specifying = £)f =1 -^H x , x ') with 
Pi(x,x') acting as P(x,x') = \x)(x'\ on the one-particle space 



p A (x,x') = trace 

q=l 



Y,^(x,x')U q R A (q) 



»=i 



(7.10) 



We want to show that, as for pi on g, off-diagonal long range order is carried by infinitely long 
cycles. For this consider the case of a positive pair potential and let H A (q) = ||p«|| 2 /2M 
be the kinetic energy of q particles. Then by dropping the interactions of the particles in 
the g-cycle with those belonging to other cycles qj, j = 1, . . . , r, one has the inequalities 

H A (q, q 1} . . . , q r ) < H°(q) + H A (q x , ...,q r ) 

trace ^...^ [U qi U q2 ■ ■ ■ Uqr e^ H ^^-^] < 

e-««)trace n m^...^ qr [U Ql U Q2 ■ ■ ■ u^-^^,..^ {7 U) 



The above traces involve off-diagonal configurational matrix elements because of the presence 
of the permutation operators U qj but the inequality is easily seen to hold using the Feynman- 
Kac path integral representation of these traces. When inserted in ()7.9|) this shows that each 
term of the sum (j7.10|) is bounded by 







q 


trace 




/ x | e -/3^°| x 



y^^Pijx, x')U q e 

i=l 

M \ d/2 ( M||a;-:r'|| 2 \ 

^m) ex H — wpa—) ■ (7 - 12) 

For the last upper bound we have used that Dirichlet boundary conditions increase the kinetic 
energy and written the explicit form of the kernel of e~^ H °^ in infinite space. Therefore 
each term of the sum (j7.1())) tends to zero uniformly with respect to A as \x — x'\ —>■ oo and 
we conclude that only the tail of the series ()7.1()|) will contribute to off-diagonal long range 
order, namely 



oo 



p c = lim lim lim V trace „ m 

Q— >oo la— m'l— >oo |a|^oo A 
q=Q 



i=l 



(7.13) 



At the moment we do not have a proof of the equality of the expressions f!7.13|) and (j2.21|) 
for a Bose gas with pair potential interactions. 

In numerical calculations of the path integral representation of the Bose gas jl()| . it is more 
convenient to use periodic boundary conditions, i. e. A is a rf-dimensional torus. Cycles with 
q particles correspond to Brownian paths on the torus with typical length ~ ^fq. To such a 
path one can associate a winding number, the number of times that the path wraps around 
the torus. Then the typical length of a winding path is at least of the order ^fq ~ (Aj 1 ^, 
the linear dimension of A and hence the number of paths such that q > |A| 2//d is larger than 
the number of winding paths. For a particle configuration , let N w be the total number of 
particles found in winding paths. Since qp A {q) is the average number density of particles in 
g-cycles, one has for A sufficiently large 

00 00 , . 



q=Q+l q>\A\ 2 / d 



If the average number density of particles in winding paths (N W )/\A\ in is different from zero 
in the infinite volume limit, we conclude that pi on g is non zero and thus Bose condensation 
occurs (provided that pi ong 7^ is equivalent for the other possible definitions of the con- 
densate density). This criterion can be compared with that for superfluidity involving the 
fluctuations density of the winding number [HI]. Thus the statistics of winding paths will 
make it possible to estimate both the superfluid density and the Bose condensate. 
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